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ABSTRACT: We investigate the effects of couplings between curvature and isocurvature perturbations before and around 
horizon-crossings during cosmological inflation. We consider a generalized two-field inflation model, in which the non- 
canonical kinetic term allows us arbitrary sound speeds of curvature and isocurvature perturbations. By using the field- 
theoretical perturbative analysis, we calculate the cross-spectrum between curvature and isocurvature perturbations and 
the corrections to curvature and isocurvature power spectra due to the presence of couplings between them. Our analysis 
confirms previous results that the cross-correlations are generated and amplified when perturbations cross the horizons. 
Moreover, we find the cross-correlation, which was previously shown to be first-order in slow-roll parameter, can be 
enhanced when the sound speed of isocurvature perturbation is much smaller than that of the curvature perturbation. This 
is because in this case the isocurvature perturbation exits its horizon much earlier than the curvature perturbation and acts 
as a nearly constant source on the curvature perturbation. 
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1. Introduction 

It is believed that the large-scale structure in our universe grows up from the primordial quantum fluctuations during a 
period of cosmological inflation (see e.g. [^] for a review). The predictions of inflation have been supported by current 
observational data [[!]]. 

The simplest model for inflation is based on the picture that a single scalar field rolls down its potential, making the 
universe inflate and also generating quantum fluctuations. However, various alternatives are investigated extensively, one 
of which is multi-field inflation models §||gg|||0|llj|0g|l|||^||||§|6[|l|@0|| 



\12\ |27| , [34|]. The goal of the studies of multi-field models is two-fold. Firstly, many inflation models based on particle 
physics or string theory usually involve many scalar fields, which can also have non-canonical kinetic terms. Secondly, it 
has been clear that any detection of primordial non-gaussianity would rule out the simplest slow-roll single field inflation 
models 1 . On the other hand, multiple field models provide us more possibilities and have been discussed extensively 
[ BQ , p3[ pl| , p2[ p5| . Thus, it is natural and important to study multi-field inflation models in details. 

However, in the context of multi-field models, except for a few specific models, even the predictions for the spectra 
of primordial perturbations are a non-trivial task. The main reason is that, there are couplings between adiabatic mode 
and entropy mode(s), even at linear level. A well-known result is that the curvature (or adiabatic) perturbation can evolve 
on super-horizon scales in multi-field inflation whereas it is conserved in single-field inflation. This is due to that the 
entropy (isocurvature) perturbation modes act as a source term in the evolution equation for the curvature perturbation. 
This phenomena was first emphasized in Jl4|]. The production of adiabatic and entropy modes for two-field models with 
a generic potential was studied in [ p"5| ] where a decomposition into instantaneous adiabatic/entropy modes was firstly 
introduced. Multi-field models with non-canonical kinetic terms have been investigated in the slow -rolling approximation 



'Meanwhile, there indeed exist various secondary effects which can also produce a detectable level of non-Gaussianity in the observed Cosmic 
Microwave Background and Large-scale Structure. 



i 



in [|11| |12J, |22|], where the adiabatic/entropy decomposition technique was also extended in [|19[ |20p w i m non-canonical 
kinetic terms. 

As has been stressed above, one of the difficulties in multi-field inflation researches is that, one cannot trace back the 
adiabatic mode to one of these scalar fields, and the entropy modes to remaining scalar fields. In general, all relevant scalar 
fields are mixed together to give one adiabatic modes and N—l entropic modes. On the other hand, this is also the reason 
why we should expect cross-correlations between adiabatic and entropy modes [^,0, B]. However, in most of the previous 
works, quantum cross-correlations between adiabatic and entropy modes before and around horizon-crossing (QaQs) r 
are expected to be small, based on the observation that the cross-correlations are of order slow-rolling parameters before 
and around horizon-crossing. This cross-correlation has been studied analytically in [25 1 as a phenomena of oscillations 



between two perturbation modes and also been investigated in details in [ 261] with canonical kinetic term and in | |27[ ] with 
non-canonical kinetic term. Numerical studies were also presented in ||28, |2^]. 

The goal of this work is to study the cross-correlations in details. The analytic treatments in [ ^B] , |2(| [l7| | are based 
on the "diagonalization" of the coupled system: time-dependent orthogonal matrices are introduced to abstract the ap- 
proximately decoupled degrees of freedom which should be quantized. Indeed, this is the standard treatment to a coupled 
system. In this note, we take a slightly alternative approach — that is, we treat the couplings between adiabatic and 
entropy modes as "two-point" interactions, and use standard field theoretical perturbative methods to calculate the cross- 
correlation and also the corrections to adiabatic/entropy spectra themselves. More precisely, we split the full quadratic- 
order action into a "free" part in which adiabatic mode and entropy mode decouple with each other and a "two-point" 
interaction part, where the two-point couplings are treated as interaction vertices. Generally speaking, this approach 
supplies us a systematic perturbative procedure to study the coulings between adiabatic and entropy modes in details, es- 
pecially in the cases where the "diagonalization" of the coupled system cannot be done easily 2 or the "time-independence" 
of the diagonolization matrices is not a good approximation. 

We consider a generalized two-field inflation model, as described in the next section. The prototype of this form of 
Lagrangian was proposed in [ pi , 35 1 and includes multi-field /c-inflation p2[ |4(i[ ] and two-field DBI model [ pi] , [32"[ [33| , 



p~jj |37| , |44| ] as special case and thus deserves detailed study. Actually, non-gaussianities in multi-field models with 
non-canonical kinetic terms have been extensively investigated in [|l||5[ ^ fu], f£|, ||[ ||[ |fj |j| ||, ||, fP| (see 
also @ mHHH, ||]). Although the main task in this note is not to study a complex multi-field model, this 
generalized Lagrangian can make our analysis of the cross-correlation in a more general background. Especially, the non- 
canonical kinetic term of our model allows us arbitrary sound speeds of adiabatic and entropy perturbations, which are 
essential for our following analysis. Our work can be viewed as generalization of the analysis in pB] , p^ , |27|] to a general 
class of two-field models with non-canonical kinetic terms and arbitrary speeds of sound for adiabatic and entropy modes, 

This paper is organized as follows. In the next section we introduce a generalized two-field inflation model, and 
describe the scalar perturbations of it. The third section is devoted to investigate the cross-correlations in detail, based on 
the field-theoretical perturbative approach. The last section is devoted to conclusion and discussion on the limitation and 
possible extension of this work. 



2. Generalized Two-field Inflation Model 

In this work, we consider a very general class of two-field inflation models with action of the form: 




(2.1) 



where X = X\ = G U X IJ and Y = XjXj with X IJ = ~\g^ d^ 1 d v (j) J , G u is the metric for the field space and 
M p = l/v87rG is the reduced Planck mass which we set to unity in the following. In two-field case, all higher order 



2 Especially, in models as considered in this paper with different adiabatic and entropy speeds of sound, i.e. c a ^ c e , the diagonalization is a non- 
trivial task. Since in this case, the system is equivalent to a coupled oscillator system with different (free-theory) energy eigenvalues (c a k ^ c e k) plus 
"time-dependent" interactions. The simultaneous digonalization of both the free-theory Hamiltonian and the time-dependent interaction is not trivial, 
and in this case the traditional effective method is the perturbation theory. Models considered in |^, ^| has the same c a and c e , in the words of 
quantum oscillators, the couple two-state system has degenerate (free) energy eigenstates, in which the diagonalization can be done easily. 



contractions among Xjj's can be expressed in terms of X and Y, e.g., 



x^x^x^xf = -^x* + x 2 y + \y 2 



etc. The model (2.1) includes multi-field fc-inflation and two-field DBI model as special cases. For example, in multi-DBI 
model the Lagrangian is P = — jj^rj — lj — V(cj) 1 ) with 

V = det (G'j - 2fX'j) 

= 1 - IfGuX 13 + Af 2 X l /xj ] - 8fX l /xjXx ] + WfX^XjXgX^ . ' ' 

This expression for determinant T> is general. In this work, we focus on two-field case, thus the last two terms exactly 
vanish, leaving us effectively V = 1 - 2fG u X IJ + 4f 2 X l /xj ] . In terms of (fuf), this is just V = 1 - 2fX + 
2f 2 (X 2 -Y). 



This form of scalar-field Lagrangian in (2. 1 ) is the most general Lagrangian for two-field models and thus deserves de- 



tailed investigations. The goal of choosing such a general Lagrangian in this note is not only because recent investigations 



on non-Gaussianities in multi-field are based on some similar Lagrangian [ pl| , |35| , |40[ |41|, |42j, |32|, g3|, ^7|, [39, 43 
but also in order to see the effects on perturbations from the structure of the theory in a wider range 3 . As we will see, the 
non-canonical kinetic term supplies us two different speeds of sound for adiabatic and entropy modes which we denote as 
c a and c e respectively, which are essential for our following analysis. 

2.1 Background Equations of Motion 

In this work, we investigate scalar perturbations around a flat FRW background, the background spacetime metric takes 
the form 

ds 2 = -dt 2 + a 2 (t)5 lj dx l dx j , (2.3) 
where a(t) is the scale-factor. The Friedmann equation and the continuity equation are 

H 2 = ^ = -(2X IJ P (IJ) -P) , 

3 3 V ' (/J> ; ' (2.4) 

p = -3H(p + P), 

where and in what follows we denote P.(ij) = g ^fu , P,{ij){kl) = ox' dj dx KL etc ' ^ or sn °rt- I n me above equations, all 
quantities are evaluated on the background. From the above two equations we can also get another convenient equation 

H = -X IJ P (IJ) . (2.5) 

The background equations of motion for the scalar fields are 

P,(U)¥ + (ZHP,(U) + P,{U)) ¥ - P,J = , (2.6) 
where Pj denotes derivative of P with respect to (j) 1 : P j = 



In this work, we investigate cosmological perturbations during an exponential inflation period. Thus, from (2.5) it is 
convenient to define a slow-roll parameter for the expansion rate 

dlnH _ H _ P {IJ) ¥4> J 
e ~ dlna H 2 2H 2 ' 1 ' ' 

In this note we do not go into details of solving the background equations of motion, but only assume that the structure of 
P(X, Y, (f) 1 ) and thus the background dynamical equations permit such an exponential expansion period. 



3 Actually the Lagrangian in < jz. l[ ) was motivated from some similar models in previous investigations. For example, in ^\ ^()| a Lagrangian of 
the form P(X,(p I ) was introduced, which described a multi-field generalization of single-field fc-inflation. In a special form P(Y,<f> IS ) with 

Y = X + b ^ 2 1 (X 2 — XjjX IJ ^j was chosen in the investigation of bispectrua in two-field models. 
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2.2 Linear Perturbations 



In this work we focus on the linear perturbations. In multi-field models, it is convenient to work in spatially-flat gauge, 
where the metric (scalar sector) is unperturbed as in ( 2.3 ), and the perturbation of the system is encoded in the perturbations 
of the scalar fields, which we denote 5(f) 1 = Q 1 for short. 

In multi-field inflation models, it is convenient to decompose perturbations into instantaneous adiabatic and entropy 
perturbations jT5[ This decomposition was firstly introduced in [|5]| in the study of two-field inflation with a generic 
potential, and was extended in [12 19|, ^0|] in two-field models with non-canonical kinetic terms. This decomposition 
technique was also generalized to non-linear perturbations [ [To| in the context of covariant non-linear formalism [Tstl- 

The "adiabatic direction" corresponds to the direction of the "background inflaton velocity", for model described in 
this work, it is 



j _ 



P, {JK )4> J 4> K 



where a is defined as 



&= \IP, {JK) 4> J 4> K : 



(2.8) 



(2.9) 



which is the generalization of the background inflaton velocity. Actually a is essentially a short notation and has nothing 
to do with any concrete field. Note that a is related to the slow-roll parameter e as a 2 = 2H 2 e. 

In this work we focus on two-field case. We introduce the entropy basis e s which is orthogonal to The orthogonal 
condition can be defined as 4 



P 



J „J — 



= S„ 



m, n 



Thus the scalar-field perturbation Q 1 can be decomposed into instantaneous adiabatic/entropy modes as: 

Q 1 = e' a Q a + e[Q s . 



(2.10) 



(2.11) 



In spatially-flat gauge, the quadratic-order action for the perturbations of the model (2.1 ) can be calculated straight- 
forwardly. After instantaneous adiabatic/entropy modes decomposition, up to total derivative terms, the second-order 
action for the scalar perturbations takes the form [BU B2[ B5[ S2J 



S 2 



dtd 3 xa 3 



^mnQ Q 



1 

'2? 



d l Q m d l Q n + E mn Q m Q n - -M mn Q m Q r 



(2.12) 



with 



K = t) 



^ 2 P,(iK)(jL)e- I a ei^e J a e r ^ , 



1 



1 S am S an + — - 1 (<S, 



3 am $crn ) 



1P(ik){jl)X kl ) e^ei , 



-Mr, 



-Mue^ei + 2Mi J e I m e J ri 



P 



(U) 



' 2P{IK)(,JL)X 



KL ) e 1 e J 



(m, n 



<J,S 



(2.13) 



where 



-Mi j 



In (2.13) we introduce 



P, 



(KI),J 



KL V NM 



P: 



\NI)(KL)P(MJ) 



P IJ+ (X MN P {MN) 



2X MN x P Qp 



(MN)(PQ) 



3H Z 



X 



KL 



-Jp-P{KI)P{LJ) 



(Pi - 2X^P {MNhI ) P {KJ) j>« + 1| (!tx"l 



^(KI)P(LJ) 



P 



2XP 



,Y 



P X + 2X (P XX + 4XP.XY + 3Py + AX 2 Py Y ) 

Px 



P x + 2XP, Y 



(2.14) 



(2.15) 



4 In specified models, other choices of orthogonal conditions are possible. The idea is to make the kinetic terms of the perturbations decouple. The 
final results for curvature/isocurvature perturbations is independent of different choices of orthogonal conditions. 



which are the propagation speeds of the adiabatic mode and entropy mode respectively. From ( |2.13| ), the kinetic term 
K mn Q m Q n has been diagonalized, as a result of adiabatic/entropy decomposition. 

For our purpose in this work, the different speeds of sound for adiabatic and entropy modes are essential for the 
investigation of cross-correlations. Actually in multi-field models, it is generic fact that c a ^ c e which was firstly point 
out apparently in [ 23 24[] in the investigation of brane inflation model. In multiple fc-inflation with Lagrangian of the 
form P(X, (j) 1 ), the adiabatic mode propagate with sound speed c s while entropy modes propagate with the speed of 
light p^ , fiQ| l. While subsequently in [^T], ^5[], it was shown that in multi-DBI models, the adiabatic mode and entropy 
modes propagate with the same speed of sound. Now it becomes clear that, in multi-field inflationary models, adiabatic 
mode and entropic modes in general propagate with different speeds of sound c a and c e , which depend on the structure 
of specific theory [ plj , [35[ ] (see also [[35} |32|, ^7| [58| |IJ ^TJ for extensive investigations on general multi-field models 
with different c a and c e ). 

It is now convenient to introduce the canonically normalized variables 



-Qa 



Qs 



a 

— i 

C e 



(2.16) 



After straightforward but tedious calculations, the quadratic action for Q a and Q s (after using conformal time r\ defined 
by dt = adrj and up to total derivative terms) takes the form: 



S^IQct, Qs 



d-qd x 



1 



where 



1 

a 



—O 2 

z ^ a 



Q?-c 2 e 



a 
a 



2 2 

a n s 



(2.17) 



z = a 



C a H 
Pss 

Pi 



(2.18) 



a 2 R 



P\ 



with 



P.s = P.iet 

Pi = Px - 



J Pi C e , 

2XPv . 



2Pi c 2 a a 2 Pl 
P.ls = (9jPi)e^V^Ce> 



2PlsP,s 

Pi 



P ss = c\Pi (ViVjP) e\e 



s ' 



(2.19) 



and R is the Ricci scalar of field space metric G/j, T>j is the covariant derivative associated to Gu (i.e. VjVjP = 
Pjj — TfjPji)- Note that these various parameters are evaluated on the background. In general, the time-dependence of 
these various parameters are complicated. In this note, in order to proceed, we introduce several slow-varying parameters: 



dln< 



Ve 



din c a 
din a 



d In c e 
din a 



din Pi 



din a 

where a is the scale-factor. 

In canonical quantization procedure, the quantum fields are decomposed as 



d In a 



dln£ 
din a 



(2.20) 



r (fe, i]) = a k u a (k, rj) + a^_ k u* (fc, 77) , 



,(k,T]) = b k u s (k,rj) +b' l _ k u* s (k,f]) , 



(2.21) 



The equations of motion for the mode functions u a and u s can be get from varying ( 2.17 ): 

{my 



2 7 2 * 

c a k - — 



elk 2 - — + a 2 n 2 j u s + H£u' a 
These two equations form a closed system for the scalar perturbations. 



z 

z'm 



(2.22) 



3. Perturbative Analysis 



As described in the Introduction, the idea in this paper is to treat the coupling between adiabatic and entropy modes as 
"two-point" interaction vertices, and to use field theoretical perturbative approaches to evaluate the cross-correlations. 

3.1 Interaction Hamiltonian 



The first line in ( |2.17| ) describes a decoupled two-field system, where the two decoupled modes can be quantized inde- 
pendently. While the second-line can be identified as the "two-point cross-interaction vertices": 



S c [Qa,Qs\ = I d v d 3 x ( -HtiQ' a Q s + — 



(3.1) 



where the dimensionless cross-coupling £ is given in (2. 18). In the operator formalism of quantization, interaction Hamil- 
tonian is needed. The Hamiltonian density which is defined by H = Tr a Q' a ~ £ can be split into two parts: H = Ho + H c , 
with 

no ^ \nl + \cl (dQ„) 2 - f z Ql + + \c 2 (dQ s ) 2 - \ - a^l - H 2 £ 2 ) Ql , 

H c = HS, 7T a Q s - ^-^QaQs , 



(3.2) 



where Ho describes decoupled system while H c describes the cross interactions. From Ho, the free-theory canonical 
momenta are (in interaction picture) are related with time-derivatives of the fields as 



dHo 



ft dHo 

Qs = — = 



thus in the interaction picture, the cross-interaction vertices can be written in terms of Q m and Q' m as 

H c = H£,Q' a Q s -QaQs ■ 



(3.3) 



(3.4) 



Since Q a and Q s are the canonical variables for quantization, the corresponding mode functions in ( 2.21 ) satisfy the 
decoupled ("free-theory") equations of motion: 



2 7 2 

c„k 



U a = , 



elk 2 



(3.5) 



a 2 fi 2 ) u s = . 



Up to the first-order in slow-varying parameters, the mode solutions with proper initial conditions are (See Appendix 
^| for details) 

u„ (V, k) = ^ (l + |) e% V=ijff£> ((1 + s a )x) , 



Us (v,k) 

with x = —c a ki] and y = —c e krj, and 



r(i + y)^)^«((i + Se) ,), 



3 1 ,„ 
v a = - + - (2e + T} e + s a ) , 

3 . / 
v s = v + 7] p — + s e v + e\v 



3 

4j> 



where v = y | — is the Hankel function of the first kind. 

The "decoupled" two-point functions for Q a and Q s are defined as 



Q s (ki,rii)Q s (k2,T]2, 



(0) 



(0) 



{2^f5 z {k l + k 2 )G kl { m , m ). 
(2ir) 3 d 2 (k 1 +k 2 )F kl ( m ,r l2 ). 



(3.6) 



(3.7) 



(3.8) 



where the supercript "(°)" means in evaluating the above expressions the coupling between adiabatic and entropy modes 
are neglected, and 



Gfe(r?i,r? 2 ) = «o-(»7i,A;)tt*(?72,k) , 



F k {m,m) = u s (rjx, k)u*(r) 2 , k) , 



(3.9) 



where u a , u s are given in (p.6|), and * denotes complex conjugate. 

In comoving gauge, the perturbation Q a is directly related to the three-dimensional curvature of the constant time 
space-like hypersurfaces. This gives the gauge-invariant quantity referred to the well-known "comoving curvature pertur- 
bation": 

K=^-Q a , (3.10) 



where & is defined in (2.9). The entropy perturbation Q s is automatically gauge-invariant by construction. In practise, it 
is also convenient to introduce a renormalized "isocurvature perturbation" defined by 5 



S=-Q s 
a 



(3.11) 



It is thus well-known result that the power spectra for curvature perturbation and isocurvature perturbation around their 
respective sound horizon-crossings are (up to the first-order in slow-varying parameters) 



= (1 - 2e - 2s a ) F Va ((1 + s a ) x) 



n 

Vs\y*) = Ps(l~2e- 2s e ) F V3 ((1 + s e ) y)\ 



c a k/aH=l > 

C<s k/aH=l ' 
2 



(3.12) 



respectively, where the various parameters are defined in (2.2C) and Viz = (J^-) t^— anc ^ ^ s = (f") 2«r am asymp- 
totic values for the power spectra on superhorizon scales, and 



F v {x) 



Hi 1] (x) 



(3.13) 



In ( 3 . 1 2 ), quantities on the right-hand-side of the equations are evaluated at the time of adiabatic or entropy sound horizon- 
crossings, i.e. c a k/aH = 1 or c e k/aH — 1, respectively. In general since c a ^ c e , adiabatic and entropy modes cross 
their respective sound horizons at different times. For later convenience, we introduce x = —c a fo] and y = —c e krj, and in 



( 3.12 ), and are their respective values around sound horizon-crossings, up to the first-order in slow-roll parameters 
which read 

x* = —c a kri\ Cak / aH=1 ~ 1 + e\ 
1 



\c a k/aH=l ) 



(3.14) 



\c e k/aH=l 1 



\c a k/aH=l 
V* = ~ C ekv\ Ce k/aH=l 

where e is the slow-roll parameter define in (|2~^). 

In the above derivation, fij, / H 2 and thus 3 — 2v s are not supposed to be small. While in the following discussion, we 
assume that 3 — 2v s is of order ~ 0(e). Actually if the perturbation mode has an effective mass comparable to the Hubble 
scale H, its quantum fluctuations on wavelengths larger than the effective Compton wavelength would be suppressed, 
then the system can be described effectively by a single-field. 

3.2 Cross-correlations 

Now we are at the point to evaluate the cross-power spectrum. As has been stressed, the idea is treat the cross-coupling 
terms as interaction vertices. From (^4), there are two types of two-point cross-interaction vertices, as depicted in fig.|l|. 
In cosmological context, perturbative calculations of the correlation functions are based on the "in-in" formalism (see 
Appendix ^| for a brief review). The leading-order cross-correlation involves one cross-interaction vertex (see fig.[j]): 

drf (Qcr (v,ki)Qs (r), fc) H c (rf) 



Qa {riM)Qs fo.fca)) =-2Ki 
= -2^i 



dri' 



m {rf) ^G kl (v, rf) F kl ( V , rf) - G kl (v, rf) F kl (77, rf) 



(3.15) 



5 There is an ambiguity in normalizing the entrop y pert urbation. Traditionally one can choose the normalization condition to ensure that "P-k* = Ts* 
when modes cross the Hubble horizon. Our choice (3.11) corresponds to Viz* I Vs* — c e /c a , i.e the ratio of the speeds of sound of isocurvature and 
curvature perturbations. Here Viz* = r P-jz(x f ) and Vs* = T^siy*)' see (3.12). 
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Qa k\ k 2 Qs Q<r ki k 2 Q s 
• 'x « • • »■ • X « • 

n i' n n i' i 

Figure 1: Diagrammatic representations for cross-correlation vertices. Black line denotes adiabatic mode and dashed line 
denotes entropy mode. A "®" denotes cross-interaction vertex. There are two type of cross-interaction vertices. A black 
dot denotes temporal derivative. 



where G, F are defined in (3.8 ). From now on, we take the massless limit (i.e. v<r = v s = |) for the decoupled two-point 
Green's functions G and F defined in (3.6)-(3.9) and also treat H etc. as constant in evaluating the cross-correlation and 
the corrections to adiabatic/entropy spectra. This is because not only that the exact Green's functions are rather difficult 
to deal with analytically, but also that the differences between using the massless Green's function and using the exact 
Green's functions for our calculations are higher-order in slow-roll parameters. More precisely, we thus treat £ on the 
same footing as the other slow-varying parameters, and identifies terms such as e£ as higher-order quantities which can 
be neglected. At the end of our calculations, the time-dependence of various parameters such as H should be taken into 
account in order to get the correct tilts of the spectra. 

The cross-power spectrum for Q a and Q s can be defined as 



a 2 (?7) 



(77, fci) Q s (n, fc 2 ) ) = (2n) 3 S 3 (fei + fe 2 ) C as (77, fci) , 



From (3.15) and (3.6)-(3.9), after a straightforward calculation, C as can be written in the form 



C as (rj,k) 



H 2 



2c a k^ 



tr c (x,\). 



with x = —c a kr] again and 



A 



(3.16) 



(3.17) 



(3.18) 



is the ratio of sound speeds of entropy and adiabatic modes, and 

_ (1 + A) 



r c (x, A) = ' J~ ' \2 -x[2 sin ((1 + A) x) Ci ((1 + A) x) + cos ((1 + A) x) (tt - 2Si ((1 + A) x))} 
2 A v 



2A 2 

1 - Ax 2 



^- [2 cos ((1 + A) x) Ci ((1 + A) x) - sin ((1 + A) x) (tt - 2Si ((1 + A) x))\ } 

1 + A j 



(3.19) 



Here we keep the ^-dependence in the expression for T c explicitly. It is not only because that there are inflation models 
where the perturbation spectra evolves quickly even after horizon crossing and never reaches the asymptotic values on 
superhorizon scales (x — > 0), but also allows us a more precise estimate of the spectra around the horizon crossing. It is 
also interesting to note that the cross-power spectrum depends explicitly on the ratio of the sound speeds for adiabatic and 
entropy modes A = c e /c a . Especially, the factor T c depends only on the ratio A, while not on c a or c e themselves. The 
dimensionless cross-power spectrum between 1Z and S is given by 6 



_*1 

2tt 2 



C !TS (r 1 ,k)=Pn^r c (x,\). 



(3.20) 



where Viz = (■¥-) is asymptotic value for the dimensionless power spectrum for the comoving curvature perturba- 
tion on superhorizon scales as before. 

(p7|)-(p^)and(^20|) are one of the main results in this note. The key point is that, at leading-order the cross-power 
spectrum is of order ~ however its amplitude is determined by the factor T c (x, A). The dependence of T c on x and A 
is depicted in fig.||. 

Two comments are in order: 



6 Our result can be compared with (e.g.) Eq.(79) in [^|. Actually one may verify that when setting c a 
Appendix ^ for details. 



!2(j ) reduces to Eq.(79) in See 




Figure 2: r c as function of x (left) and as function of A evaluated around adiabatic sound horizon-crossing at x w 1 (right). 

• From the left panel in fig.||, when modes are deep inside the horizons, r c -C 1 and the cross-spectrum are indeed 
small. This confirms previous result that cross-correlations among different perturbation modes are always negligi- 
ble when modes are deep in side their respective sound-horizons [|5j |27|], since there the system reduces to a 
collection of weakly-coupled oscillators. However, as firstly pointed out in []25[], as long as the modes get closed to 
the horizon(s), the couplings and thus cross-correlations among different modes become more and more important, 
and the cross-correlations are generated when modes cross their horizons. Our analysis also confirms this result. 
As depicted in the left panel in fig.^, when modes get closed to the horizon, T c starts to increase, and its amplitude 
is determined by A = c e /c a , i.e. the ratio of the sound speeds of isocurvature and curvature perturbations. 

• It is more interesting to note from the right panel in fig.|]that, the value of T c and thus the cross-correlation (around 
adiabatic sound horizon-crossing) can be enhanced by small A, i.e for models with c e -C c a . This phenomenon can 
be understood intuitively. As we know, the smaller the sound speed is, the earlier the corresponding perturbation 
mode crosses its sound horizon. The small c e /c a ratio implies that the isocurvature perturbation exits its sound 
horizon much earlier before the curvature perturbation exits its sound horizon. Thus in the process when curvature 
perturbation gets closed to the adiabatic horizon, the isocurvature perturbation is already well outside its entropic 
horizon and behaves as a nearly constant (rather than highly oscillating) background, which acts as a nearly constant 
source on the curvature perturbation. More precisely, the smaller A is, the longer that the isocurvature perturbation 
behaves as a source on the curvature perturbation, and the more significant this "accumulative" effect is. This fact 
causes the amplifications of both cross-spectrum between curvature and isocurvature modes and also the corrections 
to the spectra of curvature and isocurvature perturbations by small A, as we will show in the following subsection. 

3.3 Corrections to Spectra of Curvature and Isocurvature Perturbations 

Now we would like to investigate the leading-order corrections to the power spectra of curvature and isocurvature pertur- 
bations, due to the presence of the cross-interactions. It is interesting to note that the leading-order corrections to P a and 
P B from the cross-interaction vertices involve two cross-interaction vertices and thus are of ~ £ 2 , as depicted in figj^. 
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Figure 3: Diagrammatic representations of the leading-order corrections to P a and P s . Recall that there are 
two types of cross-interaction vertices, and thus there are actually four different contributions, which we do 
not show here explicitly. 



The leading-order correction to the adiabatic power spectrum can be denoted as 

(Qct (V, fci) Q<y (V, fc 2 )) (2) = (2tt) 3 S 3 (kx + k 2 ) P?) (77, k x ) 



(3.21) 



here a superscript "( 2 )" denotes that the contribution involves two cross-interaction vertices, and 
P^ (77, kx) = - 4£ 2 x / d m / d m F kl ( m , 772) 



<J —OO J — OO 

rf -G fcl (77,771) -j^-Gfej (»?,J/2) + — -y-G tl {rj,rjx)G kl (mm) 



mm dr)i 



dm 



77x7/1 d77! 



+ ~^—G kl (77, 771) -^-G fel (77, 7/ 2 ) + ~^G kl (77, 771) G fel (77, 772) 

77^772 0?772 jtf?^ 



2C 



drji / d^Ffc! (771,772) 
) J —00 

Id . <i „ , . Id 
-1— Gk, (771,77 — G fcl 77,772) H 2 — 

771772 drji dr\ 2 mm d m 



Gkt (771,77) G fcl (77,772) 



+ ^^ Gkl d^^ 1 ( ~ 7? ' ??2 - ) + ~^tf Gkl < - ?7l ' 7? ' ) Gfcl 



(3.22) 



Here G fe (771,772) = aij ^ a{m) Gk(m,m) and F k ( m ,f] 2 ) = a(v ^ a{r]2) Fk{m,m)- After a straightforward calculation, 

(2) 

Pct can be written in the following form 



Pj 2 >=P ff £ 2 r ff (*, A) 



(3.23) 



where P CT = and 



r CT (a;, A) =^rg- {a; [-16(A + l)Ci (2a;) cos (2a;) + 8(A + 1) (tt - 2Si (2a;)) sin (2x) + 16A + tt 2 + 24] 

+ [16(A 4- l)Ci (2x) (x cos (2a;) - 2 sin (2a;)) + x (tt 2 - 8(A + 1) (tt - 2Si (2a;)) sin (2a;)) 

-16(A + 1) (tt - 2Si (2a;)) cos (2a;) - 4tt cos (x(A + 1))] 
+ 4 (a; 3 + x) Ci 2 ((A + 1» - 8 (a; 2 + l) Ci ((A + 1» sin (x(A + 1)) + 4 (a; 3 + a;) Si 2 ((A + 1» 



+4 (a; 2 + l) Si((A + l).x) [2 cos (a;(A + 1)) — ttx] — 4tt cos (a;(A + 1))} + ~^9<y ( x t ^) 



with 



+ °° fzA + i) e^ -1 ) 2 

rfz 1 % [Ei {-iz (1 + A)) + l7 r] 



<7 CT (a;, A) = 3fj|(x + 7:) 2 e l2a 
Similarly, the leading-order correction to entropy spectrum is defined as 

(Q s (77, fei) Q s (77, k 2 )) (2) = (2tt) 3 <5 3 (fei + hi) Pj 2) (77, fci) , 

and we have 



pW=P s er s (y,A) 



(3.24) 
(3.25) 

(3.26) 
(3.27) 



where P s = and y = -c e kr], and 
1 



T * (»> A ) =^TF- {y [ 4 + !) Ci2 (C 1 + A_1 ) ^ + 4 + 1) Si2 (I 1 + v) ~ ^ (tf 2 + Si ((l + A_1 ) 2/) 



8A 3 y 



+7r 2 7/ 2 + 8A 2 + tt 2 ] - 8 (y 2 + 1) ACi ((l + A" 1 ) y) sin (y (A" 1 + l)) 
-4 (t; 2 + 1) A [tt - 2Si ((1 + A- 1 ) y)\ cos (y (A" 1 + !))} + ^g s (y, A) , 



with 



+00 



g s {y,\) = ^\{y + i) 2 e a y / dze'^ 1 - 1 )'—^ [Ei (-iz (l + A" 1 )) + in] 



(3.28) 
(3.29) 



The x, y and A-dependence of T a and r s are depicted in figQand fig.^]. 
Several comments are in order: 
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gure 4: r CT as function of x (left) and as function of A evaluated around the adiabatic sound horizon-crossing at x « 1 (right). 
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igure 5: Y s as function of y (left) and as function of A evaluated around the entropy sound horizon-crossing at y ~ 1 (right). 



The left figures in fig.[|and fig.|]show T a or T s as functions of x or y respectively, for different values of A. We take 
r ff as example. Since x = —c a kr\, it implies that for mode with fixed k, when deep in the sound horizon (x ^> 1), 
corrections to the power spectrum from the cross-interactions are small. Again, this verifies previous argument 
that adiabatic and entropy perturbations can be treated as decoupled when they are deep inside the horizon, while 
when modes approach the horizon, i.e. x w 1, the modification to the power spectrum starts to increase [g5J. The 
conclusion is the same for T s . 



The most interesting point is that, as for the cross-correlation between adiabatic and entropy modes, the strength 
of the corrections to adiabatic/entropy spectra due to the presence of cross-couplings are also determined by the 
parameter A = c e /c a , i.e. the ratio of the sound speeds for entropy and adiabatic modes. It has been known that 
for c a = c e i.e. A = 1, the cross-correlation and also the corrections to the "decoupled" spectra are proportional to 
the cross-coupling and thus are expected to be small [^5 26 27|] . However, from the right figures in fig.[| and fig^, 
this corrections can be enhanced by small A, that is for models with c e <^ c a . Especially, this enhancement is most 
significant for the adiabatic mode (fig.Q). As explained before, when c e <C c a the entropy mode cross the horizon 
much earlier than the adiabatic mode, and thus act as a nearly constant (rather than highly oscillating) source on the 
evolution of adiabatic mode. Thus the smaller A is, the longer that the entropy mode behaves as a source, and the 
more significant this accumulative effect is. 

Inversely, this enhancement is not significant for the entropy mode. Actually from the right panel in fig.||, when 
A > 1 i.e. c e > c a , the correction to the power spectrum of entropy mode is suppressed when A goes large. 
Intuitively, this is because that, as is well-known, on super-horizon scales entropy modes can act as sources for the 
evolution of adiabatic mode, while inversely adiabatic mode can never act as a source for entropy modes. Thus, the 
larger A is, the earlier the adiabatic mode exits its horizon, and the less it affects the evolution of entropy modes. As 
an extremal case, one can verify that T s — > when A — > oo, which implies that in this case the entropy mode is not 
affected by adiabatic mode and evolves freely. Moreover from the right panel in fig.Bl when A — > 0, T s approaches 
a constant value rather than blowing up, since in this case the adiabatic mode which is highly oscillating affects the 
entropy mode with a nearly constant strength. 



What we are eventually interested in are the spectra of the curvature and isocurvature perturbations, which are defined 



in ( |3.10| ) and ( |3.11| ). After using ( |3.12| ), the power spectrum for the curvature perturbation, including the leading-order 
corrections from the cross-interactions, and also including the corrections to both Green's function and Hubble parameter 
in first-order slow-roll parameters, takes the form 



■R 



1 - 2e - 2s a + [v a - - / ((1 + s a ) x) + s a 



2x 2 



(1 



-r CT (x, a) 



where 
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d 



2(1 + x 2 ) \dv 



(x) 



(3.30) 



(3.31) 



In deriving the above expression, (C.4) and (|C.4|) are used. Similarly, the power spectrum for isocurvature perturbation is 



V s {y)^V s (l + y 2 



1 - 2e - 2s e + [u s - - / ((1 + s e ) y) + s 



2y 2 



i + y 2 (i + y 2 ) 



r, (y, A) 



(3.32) 



(3.20), (3.3C) and (3.32) are the main results in this note, and can be viewed as generalizations of previous results, e.g. 
Eq.(78)-(80) in [^. 

3.3.1 Deep Inside the Horizon 

In [0], an oscillating mechanism was introduced to study the cross-correlations between perturbations, where it was 
found that when deep inside the Hubble horizon different modes evolve independently and can be considered as good 
mass eigenstates, thus the cross-correlations are indeed small. Intuitively, when deep inside the horizon, the system 
become weakly-coupled oscillators in Minkowski background in which the couplings among them are assumed to be 
small (of order slow-roll parameters). In fact, as an explicit confirmation, one can show that 



r e (x,\) 



2»1 



2 + A 
(T+A) 2 



(3.33) 



which is independent of x. Similarly, it can be verified that g,j yS (x, A) approach constant values (independent of x) when 

x > 1, 

g a {x,X) -> — — — -r, 9s(yA) 



2(1 + A) ' 



2(1 + A) ' 



thus 



F<j{x, A) 

r s (y,A) 



1 



2(1 + A) 2 ' 
1 + 2A 



(3.34) 



2(1 + A) 2 



From (3.33) and (3.34), it is obvious that 



maxTi - 0(1), 



(3.35) 



thus we can conclude that when deep inside the horizon, the cross-correlation between curvature and isocurvature pertur- 
bation is always smaller than O (£), and the corrections to curvature/isocurvature perturbation due to this cross-interactions 
are always smaller than 0(£ 2 ). This confirms previous investigation that the couplings between adiabatic/entropy pertur- 
bations can be neglected when modes are sub-Hubble. 

However, as was firstly pointed out in j25] ] and was also analyzed in [[l6| ^7|], the most important lesson we get is 
that, the cross-interactions between adiabatic and entropy modes at linear level 7 , which are negligible when modes are 
deep inside the horizons, are generated and amplified when modes cross their sound horizons. Our analysis in this work 
also confirms this fact. 



7 When investigating non-Gaussian features, higher-order interactions among different field modes are important, e.g. in multi-field models or in 
inflaton-curvaton models. Loop corrections also involve cross-interactions, see e.g. n59(l 



3.3.2 Around Horizon-crossings 



As has been stressed before, ( |3.20[ ), ( |3.30| ) and ( ft.32| ) are the general expressions for power spectra for curvature and 
isocurvature perturbations and their cross power. Thus in general ( 3.20| ), ( 3.30| ) and ( 3.32 ) are needed to evaluate more 
precisely the amplitude of the powers and the spectral tilts around horizon crossing (x* ~ 1 + e a and ~ 1 + e e ). In this 
subsection, for our purpose to get a glance of the effects of the cross-correlations, we use superhorizon asymptotic limits 
of (3.2C), ( 3.3C| ) and (3.32) to evaluate various quantities around horizon-crossing 8 . 
The cross-power spectrum around adiabatic horizon-crossing is approximately 9 

1 + A - ln(l + A) - 7 
....... - 2 , 



Cnsls^-i - ^C(A) with C(A) = 

aH ' a.H ~ L 

where 7 ps 0.5772 is the Euler-Mascheroni constant. The spectral index is 

-2e - i] + 775 + Cs a 



(3.36) 



d\n Cks 



din A; 



(c + 1 ) 



-=1 



(3.37) 



with C = 1 — 1/((A + 1) C(A)). Note that in (3.36) and (3.37) all quantities are evaluated at the time of adiabatic sound 
horizon-crossing, i.e. c a k/aH = 1. From (3.36) it is explicit that up to first-order in the cross-interaction coupling £, the 
cross-power spectrum is <~ C(A)£, it can be enhanced by small A due to the factor C(A), which scales as 1/A 2 for small 
A. Note that for A = 1 i.e. for c a = c e , C(A) reduces to the familiar value C(l) = 2 — In 2 — 7 w 0.7296. Thus our 
calculations can be viewed as generalization of previous results in j2^ , p7| ]. 

Similarly, power spectra for the curvature and isocurvature perturbations around their respective horizon-crossings 

are 

Vn\^L =1 =s 7Ml-2e-2s a ) [l + (3 - 1v a ) (7 - 2 + ln2) + £a„{\)\ \ Cak=1 , 

a.H ^ aH (3 38) 

Ps\<g =1 - V s (l-2e-2s e ) [l + (3-2^)( 7 -2 + ln2)+e 2 a s (A)]| £ ^ =1 , 

where a CTiS (A) are functions of A only. In getting ( 3.38| ), we have neglected all terms proportional to (such as) e£ 2 etc., i.e. 
we only keep the leading-order contributions from £ 2 . For general A, it is difficult to abstract analytical expressions for 
a ffjS (A). Their numerical results are depicted in fig.|[ Form which it is explicit that the corrections to both curvature and 
isocurvature power spectra can be enhanced by small A. 
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Figure 6: Numerical evaluation for a a (A) and a 3 (A) as functions of A. 

In practice, it is convenient to introduce a dimensionless correlation angle A to characterize the strength of the cross- 
correlation: 

cosA=^£Sf=. (3.39) 



From (3.20), (3.3C) and (3.32), and using the fact that around the adiabatic sound horizon-crossing V-jz/Vs ~ A, we get 
the correlation angle around adiabatic sound horizon-crossing 



cosA|c< l fc =1 £VXT c (x, A) 



(3.40) 



This is the traditional treatm ent in t he lite ratur e Ho wever, the spectral indices of the spectra in multi-field models in general should be evaluated 
by using the general expressions (3.20), (3.30) and (3.32), rather than naively by using the curvature perturbation as in single-field models. See e.g. the 
numerical discussions in [ p7| . 

9 ( 3.36l can be compared with (e.g.) Eq.(37) in [po|. If we set A = I and § = —2r) as , (3.36) reduces to Eq.(37) in [E3]. 



1 1 



In getting (3.40) we take the approximation under the assumption that the power spectra of curvature/isocurvature pertur- 



bations are still dominated by their "decoupled" values V-j^ an d T~'s > ^ i- e we assume that the corrections to the spectra 
would not exceed their respective decoupled values (they are indeed corrections). A more precise evaluation can be done 
by using (3.36) and (3.38) and the corresponding correlation angle around adiabatic horizon-crossing is depicted in fig.^|. 
It immediately follows that, for A = 1 (the case for canonical kinetic terms or multi-DBI case with c a = c e = c s ) or 




Figure 7: The correlation angle cos A» around adiabatic sound horizon-crossing as function of A. Here £ is chosen as 0.1 



larger (a concrete example is the case for multi-field /c-inflation [ |22| , |40j, where c a = c s <C 1 and c e = 1), the correlation 
angle is small and the cross-correlation between curvature and isocurvature perturbations are suppressed. However, as 
has been stressed before, for models with c e < c a or even c e <€. c a , the cross-correlations will be highly enhanced by 
small A. Thus, for models with £ much larger than the slow-roll parameters which are of order 1CP 2 or with c e -C c a , the 
curvature and isocurvature perturbations are highly correlated when exiting the adiabatic sound horizon. In this case, the 
"decoupled" power spectra and are not good approximations, and the cross-correlations between curvature and 
isocurvature perturbations much be taken into account. 



4. Conclusion and Discussion 

In this paper we investigate the effects of cross-interactions between curvature perturbation and isocurvature perturbations 
on cross-power spectrum and the power spectra of curvature/isocurvature perturbations themselves, before and around 
horizon crossing. Previous investigations of the cross-correlations are based on diagonalizting the coupled equations 
of motion |27j]. However, one can verify that, for models with different adiabatic/entropy speeds of sound, the 

"diagonalization" cannot be done easily and thus the treatments in pBj , |26, 27 1 cannot be viewed as good approximation 



in this case. Thus, in this work, the cross-couplings are taken as two-point interaction vertices, and a field-theoretical 
perturbative approach is taken to evaluate the cross-power spectrum etc 10 . 



The main results in this work are summarized in (3.20), (3.3C) and (3.32). Our analysis confirms previous conclu 



sion that the cross-correlations, which can be safely neglected when modes are deep inside the horizons, are generated 
when modes cross their sound horizons p5|]. Moreover, the most interesting phenomenon get in this work is that the 
cross-correlation (and also the corrections to power spectra of curvature/isocurvature perturbations) can be enhanced by 
small c e /c a ratio, where c a and c e are the sound speeds of curvature and isocurvature perturbations respectively. As has 
been stressed before, this happens since in models with c e /c a <C 1, in the process curvature perturbation getting closed 
to the adiabatic horizon, the isocurvature perturbation is already well outside its entropic horizon and behaves as a nearly 
constant (rather than highly oscillating) background, which acts as a nearly constant source on the curvature perturba- 
tion. This fact causes the amplifications of both cross-spectrum between curvature and isocurvature modes and also the 
corrections to the spectra of curvature and isocurvature perturbations. 

To end this note, we would like to comment some limitations and also possible extensions of the investigation in this 
note: 



10 This perturbative method is standard, which has also been introduced independentl y b y Chen and Wang very recently in evaluating the "transfer" 
between adiabatic and entropy modes as well as their implications on non-Gaussianities[r>lh. 
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• The approach in evaluating the cross-correlation and the corrections to (decoupled) power spectra of curvature 
and isocurvature perturbation is perturbative. The coupling £ is assumed to be small. This is indeed the case for 
canonical kinetic terms, but in models (e.g.) considered in this note, £ may be not small enough for the perturbative 
approach to be valid. 

• In general the cross-couplings between adiabatic and entropy modes are complicated. In this note we make the 



assumption that the cross-coupling £ in (3.2) is slow-varying in time. This is a simplification for abstracting the 



basic properties of the cross-correlations, but definitely a more detailed analysis of the cross-interactions is needed. 

Moreover, in this note we only focus on the cross-correlations around the horizon-crossing. In order to get the 
resulting primordial power spectra on large scales, and to compare the predictions of a multi-inflation model with 
observations, one must then solve the coupled system described by the full equations of motion. In some particular 
case and within the slow-roll approximation, one can arrive at an analytical expression for the spectra on large 
scales, e.g. through the "transfer matrix" method [ph. In general, however, a numerical approach is needed. 



As has been stressed before, ( 3.20 ), ( 3.3C ) and ( J3.32| ) are the general expressions for power spectra for curvature 



and isocurvature perturbations and their cross power. In general, in multi-field models it may not be permitted to 
use the later time (superhorizon limit) asymptotic expressions to evaluate the powers around horizon-crossing, since 
there may be inflation scenarios where the perturbation evolve quickly after horizon crossing and never reach the 
asymptotic values. Thus in general ( ^20| ), ( fOp) ) and ( |3~32| ) are needed to evaluate more precisely the amplitude of 



the powers and also the spectral indices around horizon-crossings = 1 and |# = 1 respectively). 



One possible application of the formalism and result in this work is that, as firstly pointed out in pCp, one should 
expect the "transfer" of non-Gaussianities from entropy perturbations to the adiabatic perturbation (see fig.^j), if the 
cross-correlations between adiabatic mode and entropy modes are larger than ~ 0(e) and thus cannot be neglected. 
Especially, there are inflation scenarios in which the non-linearities in adiabatic mode itself are small, however 
the possible large non-Gaussianities in entropy modes could transfer to the non-Gaussianities in adiabatic mode 
through the cross-interactions. This would bring new features to non-Gaussianities such as new shapes of momenta 



Figure 8: Diagrammatic representations of exchanging non-Gaussianities from 
entropy mode to non-Gaussianities of adiabatic mode. 

configurations [30, 6l|]. 
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A. "In-in" Formalism 

The "in-in formalism" (also dubbed as "Schwinger-Keldysh formalism", or "Closed-time path formalism") [^5j [46[ 
(see also [Q, f|9|] for a nice review) is a perturbative approach for solving the evolution of expectation values over a finite 
time interval. It is therefore ideally suited not only to backgrounds which do not admit an S-matrix description, such as 
inflationary backgrounds. 



In the calculation of S-matrix in particle physics, the goal is to determine the amplitude for a state in the far past \ijj) 
to become some state \tp') in the far future, 



= (^'(+oo)M-oo)) 



Here, conditions are imposed on the fields at both very early and very late times. This can be done because that in 
Minkowski spacetime, states are assumed to be non-interacting at far past and at far future, and thus are usually taken 
to be the free vacuum, i.e., the vacuum of the free Hamiltonian Hq- The free vacuum are assumed to be in "one-to- 
one" correspondence with the true vacuum of the whole interacting theory, as we adiabatically turn on and turn off the 
interactions between t = — oo and t = +00. 

While the physical situation we are considering here is quite different. Instead of specifying the asymptotic conditions 
both in the far past and far future, we develop a given state forward in time from a specified initial time, which can be 
chosen as the beginning of inflation. In the cosmological context, the initial state is usually chosen as free vacuum, such as 
Bunch-Davis vacuum, since at very early times when perturbation modes are deep inside the Hubble horizon, according 
to the equivalence principle, the interaction-picture fields should have the same form as in Minkowski spacetime. 

The Hamiltonian can be split into a free part and an interacting part: H = Hq + Hi. The time-evolution operator in 
the interacting picture is well-known 



U{ V2 , m ) = Texp (-if™ dt'H a (r)')^ 



(A.l) 



where subscript "I" denotes interaction-picture quantities, T is the time-ordering operator. Our present goal is to relate the 
interacting vacuum at arbitrary time |Qi(i)) to the free vacuum |0i) (e.g., Bunch-Davis vacuum). The trick is standard. 
First we may expand |Oi(??)) in terms of eigenstates of free Hamiltonian Hq, \ili(rj)) = J2 n \ n i) ("il^M 7 ?))' tnen we 



evolve |Oi(77)) by using (A.l) 



PM) = c%2,m)|fii(»7i)) = |0i) (OilOj) + $>+ 



iE n (ji2—rii) 



ni) (ni|fii(rii)) 



(A.2) 



n>l 



From (A.2), we immediately see that, if we choose 772 = — oo(l — ie), all excited states in (A.2) are suppressed. Thus we 
relate interacting vacuum at T] = — oo(l — ie) to the free vacuum |0) as 



|ni(-oo(l-ie))) = |0i) (Oi|fii) 
Thus, the interacting vacuum at arbitrary time 77 is given by 

|VAC, in) = \tk{rj)) = Ufa -oo(l - ie))|fii(-oo(l - ie))) 
= Tcxp ( -i P dr/Hntf)) ^ (0i|O t ) . 

y J-oo(l-ie) J 

The expectation value of operator 0(rj) at arbitrary time -q is evaluated as 
_ (VAC,in|0(r?)|VAC,in) 



(A.3) 



(A.4) 



(VAC,in|VAC,in) 



= 0, 



TVxpj / / ' dr/Hutf)] 01(77) Tcxp \-i f drfHairf) 

oo(l+ie) I \ J—oo(l—ie) , 



(A.5) 



where T is the anti-time-ordering operator. 
For simplicity, we denote 



-oo(l — ie) — —00^ 



-oo(l 



ie = —00 



(A.6) 



since, e.g., — 00 + has a positive imaginary part. Now let us focus on the time-order in (A.5). In standard S-matrix 



calculations, operators between (0| and |0) are automatically time-ordered. While in ( A^), from right to left, time starts 
from infinite past, or — 00 + precisely, to some arbitrary time 77 when the expectation value is evaluated, then back to — oo~ 



1 1 



again. This time-contour forms a closed-time path, so "in-in" formalism is sometimes called "closed-time path" (CTP) 
formalism. 

The starting point of perturbation theory is the free theory two-point correlation functions. In canonical quantization 
procedure, we write a scalar field as 

<t>k{v) = u(k,r])a k + u*(k, r/)a_ fe , (A.7) 

where u(k,rj) is the mode function for 4>k(v) ( m practice, Uk(rf) and u%{rj) are two linear-independent solutions of 
equation of motion for <pk(v)> which are Wroskian normalized and satisfy some initial or asymptotic conditions ). 
The free two-point function takes the form 



<O|0 fcl (j7i)0 fc2 (Tft)|O> = (2Tr) 3 6 3 (k 1 +k 2 )G kl ( m , V2 ), 



(A.8) 



with 



G fcl (7/1,7/2) = u kl {r)i)u* ki (rt2 



In this work, we take (A.8) and ( |A.9[ ) as the starting point. 
Now Taylor expansion of ( |A.5[ ) gives 

• Oth-order 



(0) 

0(7/)) = (0i|Oi(r?)|0i) 



1 st-order (one interaction vertex) 



(1) 

0{r))) =2 Re 



OO + 



dr)' (Oi 0i(r?) H a {rf) 



(A.9) 



(A. 10) 



(A.ll) 



• 2nd-order (two interaction vertices) 



(2) 



0(7?)) =-2 Re 



00+ 
n 



dr]' 
dr{ 



dr," (Oi 



00+ 
n 



dr]" (0 



6^) H^v') H a ( V ") 0j 



(A. 12) 



B. Mathematics 

In this work, we frequently account exponential/sine/cosine-integral functions, their definitions are 



Ei(z) = I —dt 

1 sin(f ) , 
Si(z) = / —j^dt, 



In evaluating various integrals, the following properties are frequently used: 
• For x > 0, 



, cos(i) , 
Ci(z) = - / — ^ dt . 



Ei 



(-ix) = Ci{x) - i (Si{x) + |) , 



where Ci(x) and Si(a;) take real values. 

• Ei(— ioo) = — in, Ei(+ioo) = in, Ci(+oo) = 0, Ci(— 00) = in and Si(±oo) = 

• When x — > +00 

e~ ix ( 1 

Ei(— ix) + in = i h 0(x 

x V x 



(B.l) 



(B.2) 



I "7 



C. Comparison with Previous Results with c a = c e 



In order to compare calculations in this note to previous results, we consider a special case c a = c e , or A = 1. Our 
goal is to show that the field-theoretical perturbative approach taken in this note is essentially equivalent to previous 

|. We take T c as example. 



"diagonalization" method taken in [|25| gg, 

It is straightforward calculations to verify that 



r c (x*, 1) = - (a^ - 1) [2Ci(2ar») cos(2a;») - (ir - 281(20;*)) Bva(2x m )] 
- x [2Ci(2a:*) sin(2x*) + (tt - 2Si(2a;*)) cos(2a:*)] + 2 

dY v (a;* 



4 X * 



(CI) 



Recall that the dimensionless power spectrum for a scalar field in dS background takes the form 

V{r),k)=VF v {x) , 

where P = (H^-) is the asymptotic super-Hubble limit (i.e. — > 0) of the spectrum and x = ~c s kr\ and 



(C2) 



F v (x) = ~x 3 (x) 



(x) 



3\ 7T 



2/ 2 



5!/ v 1 



(C3) 



= 1 + x/ 



" ~ 2 J /(*) + 



with 



2(1 + x 2 ) V^J/ 



' ? (x) 2 



(C4) 



Thus in this case the factor r c (x*, A) which is get from theoretical perturbative calculations is nothing but the coefficient 
of the term which is first-order in v — I in the expansion of the exact power spectrum: 



r c (x», i) = 



(C5) 



Thus (ZA) can be compared with (e.g.) Eq.(79) in ^Tfl. This result is as expected, since in this case the original coupled 
equations can be diagonalized, and two approaches (i.e. "calculating the correlator perturbatively" and "expanding the 
exact correlator perturbatively") are essentially equivalent. 



D. Mode Solutions 



We take the adiabatic equation of motion in (3.5) as example. It is convenient to introduce new variable x = —c a kr) 
to characterize the time-evolution 11 . After straightforward algebras and keeping terms up to first-order in slow-varying 
parameters, the equation becomes: 



S a d,U 



dx 2 



1 



x ~df + ' 1 + 2Sa ~ ~ ' 2 + 3e + " V ' 



x* 



u„ = . 



(D.l) 



where various slow-varying parameters are defined in (2.20). Its solution which corresponds to the traditional BD vacuum 
is 

UtT = Cx^ 1+s ^H^((l + s a )x) , (D.2) 



'A similar treatment can be found in Appendix A in j6c|], where an alternative variable y = was chosen as the evolution variable. 



i o 



where C is the normalization constant to be determined and v a = | + | (2e + r? e + s a ). 

On the other hand, when modes are deep inside the horizon i.e. x« 1, the above equation becomes 

d 2 u a s a du 
dx 2 x dx 

and its solution which satisfies the Wronskian normalization is 



+ (1 + 2s a )u a =0. 



= a pi(l + Sa)3 

(tf(l-e)) 2 , , .... ' 



l+s a 



(D.3) 



(D.4) 



where we use a " " to denote constant values for the slow-varying quantities, e.g. H ~ Ha e etc. Using the fact that for 
x 1, the solution (p.2|) has the following asymptotic expression: 



2 e-'ie" 



7T V / T T K 



the constant in (p.2|) can be determined as 



This gives the solution in (3.6). It is similar for the entropy mode solution in (3.6). 



(D.5) 



(D.6) 
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